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Abstract 

During solar flares fast electron beams generated in the solar corona are 
non-thermal radio sources in terms of type III bursts. Sometimes they can enter into 
the interplanetary space, where they can be observed by in-situ measurements as it 
is done e.g. by the WIND spacecraft. On the other hand, they can be the source of 
non-thermal X-ray radiation as e.g. observed by RHESSI, if they precipitate toward 
the dense chromosphere due to bremsstrahlung. Since these energetic electrons are 
generated in the corona and observed at another site, the study of transport of such 
electrons in the corona and interplanetary space is of special interest. The transport 
of electrons is influenced by the global magnetic and electric field as well as local 
Coulomb collisions with the particles in the background plasma. 



1 Introduction 

In the solar corona energetic electrons are released e.g. during solar flares and travel 
along magnetic field lines either toward or away from the Sun (Figure 1). While they 
propagate through the plasma background they excite Langmuir waves via beam-plasma 
instability (see e.g. Melrose, (1985)). Partly the energy of Langmuir waves is converted 
into electromagnetic radiation with a frequency close to the electron plasma frequency / 



/ = / (1) 

and/or its harmonics, where eo, e, m e and N e denote the electric constant, the elementary 
charge, the mass of an electron and the electron number density, respectively. 

If those electrons have sufficient energy, and the plasma background is dense enough (as 
in the case of the solar chromosphere), they can emit non-thermal X-ray radiation via 
bremsstrahlung (see e.g. Brown, (1972)), which can be detected by spacecrafts as RHESSI. 
On the other hand electrons moving toward the interplanetary medium up to 1 AU can 
be observed in-situ by spacecrafts as WIND. 
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Both methods provide the energies of the elec- 
trons either at the X-ray emission site (which is 
different from the electron acceleration site) or at 
the current position of the spacecraft in the in- 
terplanetary medium. With this information and 
the model explained in the present paper, it is 
possible to calculate the altitude of the accelera- 
tion site and the kinetic energy that the electrons 
would have had there. 

Figure 1: Left: The dynamic spectrum shows type III 
(and U) radio bursts in the solar corona, which were 
detected by the AIP spectral polarimeters. (The in- 
tensity is colour coded.) Right: The propagation of 
electrons is schematically illustrated. The fast prop- 
agating electrons excite type III radio emissions. 



In Section 2 the theoretical model to describe the electron propagation through a given 
plasma background will be introduced. Next, the plasma background itself is going to be 
explained in Section 3. Finally in Section 4 some numerical solutions will be presented 
and in Section 5 the main issues summarised. 

2 Model for Electron Transport 

The equations of motion for an electron propagating through a given plasma background 
(the latter explained in Section 3) will be briefly introduced within the current Section, 
by taking Coulomb collisions, global magnetic and electric fields into account. 

Therefore, first the assumptions are explained, which are made to obtain the equations 
of motion and next the important relations and appearing quantities. 

Assumptions 

1. Because collisions are easily to treat in the centre of mass system, all calculations 
are done in this frame of reference. This means, all quantities describe only an 
effective one-particle system consisting of the electron and the plasma background 
particles. 

2. Within this paper a spherical, heliocentric coordinate system is used consistently. 

3. The electric and the nonuniform magnetic fields are treated to be collinear and 
one dimensional. They are introduced within Section 3 and have only small spatial 
variations, i.e. the magnetic moment becomes an adiabatic constant of motion. 
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4. The acceleration by Coulomb collisions is considered to be caused under a small an- 
gle scattering approximation, due to the fact, that large angle scattering phenomena 
are rare in the cases presented within this paper. 



Important Relations, Quantities and the Final Equations of Motion 



Due to Coulomb collisions, the mean squared scattering angle 6 2 j will change, according 
to Jackson (1962), with the time t as 



d6>2 



e 2 q 2 N, 



In 



\ 



dt, 



(2) 



if the electron during its passage through the background plasma is affected by a particle 
j. Here j represents a particle of the plasma background, with charge qj and parti- 
cle number density Nj, while Vj and jij stand for the relative velocity and the effective 
mass between the background particle and the propagating electron. The other appear- 
ing quantities are the Debye length Ad = {eok^,T{N c e 2 )~ l ) Q ' b , the Boltzmann constant 
k-Q and the temperature T. b = e\qj\ (AiteQiijV 2 ) -1 characterises an impact parameter 
for the effective one particle system, which leads to a deflection of 0.5 n and allows to 
define 6 ,|| := bo cos -2 [0^] by using the pitch-angle 0j — IB, Vj . Hence the following 
root-mean-square (RMS) quantity is introduced by definition 



d0, 



(3) 



The effects caused by the electric field E and the magnetic flux density B to the pitch 
angle Qj are given, as Bai (1982) suggested, by 



dQj sm [©i] e E Vj sin[0j] dB 
dt Vj fij 2B dr 



(4) 



Because of assumption 3 stated on page 2 it is fully sufficient to investigate the scalar 
quantities only. 

If one considers more than one particle species j in the plasma background, as it is 
described in Section 3, then one has to introduce the following quantities 



d#c,total 
©total 
A^total 



d0 C ,e + d#C,p + d#c iH e2+ 

e c + e p + e Hc 2 + 



(5) 
(6) 
(7) 



where k e ,k p and /c He 2 + gi ye the divvies of electrons, protons and a-particles inside the 
plasma background i.e. j E {electron (e), proton (p), a— particle (He 2+ )}. 

The total velocity is affected both by the electric field and by Coulomb interaction namely 

as 

dftotal 



dt 



eE dw c ,totai 

COs[0 to talJ + 



/^total 



dt 



(8) 
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Again if one considers the same plasma background as above, one obtains for the Coulomb 
acceleration, which is the last term of Equation (8) the following expression 



dv, 



C, total 



dv* dvt 
+ 



+ 



dv 



He 2 + 



dt dt dt dt 

where each addend is given according to Jackson (1962) by 



d^ 
dt 



e 2 Assign 



4ne 2 



In 



>l±3 

2bl 



(9) 



(10) 



Hence, Equation 9 shows that the total Coulomb acceleration is given by the sum of 
the electron, proton and He 2+ accelerations. Under the convention that v* < for an 

electron moving toward the Sun and v* > for one propagating outward, sign v* in 
Equation (10) makes sure, that the acceleration is in fact a deceleration of the electron 
velocity. Therefore Equation (10) represents a kind of friction, which forces the electron 
to stop due to its energy loss per collision. At last, the change of the radial component r 
of the coordinate of the effective one-particle system's movement is given by 



dr 
dt 



^total COS [©total] == f total, || 



(11) 



If the initial altitude r := r[t ], initial velocity f to tai,o := ^totai^o] an d initial pitch-angle 
©totai.o := ©total [to] with the introduced previous equations are given then it is possible to 
calculate the propagation of an electron away from the acceleration site with the following 
final equations of motion: 



r[t + dt] 
^totai[£ + dt] 

©total [t + dt] 



VtotalM COS [©total [t]] dt 
dVtotal 



= UtotalM + 



dt 



dt 



t,e t 



^ W , ( d@total , ^d6 c .total 

ftotai m + i — - — i- y- 



dt 



dt 



dt 



(12) 
(13) 

(14) 



t,®total[*] 



Here V denotes the result of a binary dice, whose results —1 and 1 are equally distributed. 
This is needed to consider the random Coulomb scattering directions. Equation (14) 
explains the relation between © to tai and #c .total- These presented Equations are sensitive 
to the models of the plasma background, i.e. the plasma densities, the magnetic fields and 
the electric fields. The geometry of the electric and the magnetic field has been arranged 
by assumption 3 on page 2. 

Note that those equations allow a so-called backward calculation (dt < 0). Therefore this 
method provides a good possibility for estimating the altitude and the energies of the 
electrons both at the acceleration site (t = =: t ), if for any t f > t the altitude r[t f ], the 
velocity ^totai^f] an d the pitch-angle ©total [if] are known. 

These kind of information (as mentioned in the Section 1) can be obtained e.g. by com- 
bining RHESSI and WIND data. Then, one has to assume a standard flare model, which 
implies that the acceleration site emits the electrons in both directions (in- and outward) 
simultaneously. 
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Figure 2: Left: The model of Newkirk (1961), as given by Equation (15), is plotted for several 
different as (dashed lines) together with the introduced model by Mann et al. (1999) (solid 
line). It is clearly visible that r§ is a function of a. Right: The density model according to 
Mann et al. (1999) is plotted up to a range of 1 AU ps 215 i?sun- 

3 Plasma Background 

Within this Section the plasma background i.e. the model for the densities, the global 
magnetic and electric fields are introduced. 



3.1 Density Models 

The composition of the background plasma is assumed to be fully ionised and con- 
taining electrons, protons and a-particles: N e = k e N m 0.52iV, N p = k p N ps 0A4N, and 
-/V He 2+ = k Ke 2+N 0.04iV, where iV denotes the total particle number density, which is 
the sum of the electron N e , the proton N p and the He 2+ iV He 2+ densities. Note that these 
quantities lead to a mean atomic weight of Ji = 0.6 as given in Priest (1984). 

For the coronal density model, one can use the a-fold-Newkirk (1961) model illustrated 
in Figure 2 and given with the following equation 

AWi rk [r,a] = 4.2afc" 1 x 10 10+4 - 32 m -3 > (15) 

while i?sun ~ 6.958 x 10 s m denotes one solar radius. The a parameter is supposed to be 
about 1 in the case of the quite, about 4 for the active and about 10 for the super active 
Sun. Due to the fact, that the Newkirk (1961) model is only valid close to the Sun, an 
interplanetary density model has to be taken into account additionally. Mann et al. (1999) 
introduced such a model, which considers the stationary flow of the solar wind, by solving 
the wind equation of Parker (1981) and the equation of continuity simultaneously 

Nvr 2 = 1.2 x 10 34 Hz, (17) 
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with v c := (k B Tfi~ 1 m~ 1 ) - 5 and r c := 0.5 GM Sun v~ 2 . Here G stands for the Newtonian 
constant of gravitation, M Sun ps 1.985 x 10 30 kg represents the mass of the Sun and m p the 
mass of a proton, while v is the velocity of the solar wind. For a temperature of T = 10 6 K 
one obtains v c ~ 120 Mms -1 and r c ps 4.2 -Rsun while the numerical solution of the Equa- 
tions (16) and (17) provide the interplanetary density model of Mann et al. (1999), which 
is shown in Figure 2. 

To combine the presented models of Newkirk (1961) and Mann et al. (1999) one has to 
introduce rs as the point of intersection between both of them, then 

Mr d = { iVNcwkirk[r ' a] for li?Sun <r<rs (18) 
l ' J I N Mann [r] 5AU>r>r s ' 1 * J 

with a = 4 (and rs[a — 4] m 1.12 i?sun) yields to the density model, which is used within 
this paper (see Onel (2004,2005) for more information about this explained density model). 



3.2 Global Magnetic Field 

In the present paper the coronal magnetic field model .Bd&m of Dulk & McLean (1978) and 
the interplanetary magnetic field model -Bm&n of Marini & Neubauer (1990) by considering 
Musmann et al. (1977) as well as Parker (1958) are combined as elucidated in Onel (2004). 

This yields to a global magnetic field model which is plotted in Figure 3 and described 
by the following equations: 

3/T \ - 3 / 2 

V t / v 7 ! xl0~ 4 T for 1.02 R Snn < r < 10 i? Su n (19) 




R 



Sun 



7.988 T 2 m 4 9.6112 x 10 7 T 2 m 2 - 2 . 1A „ . . t ATT , nn . 
~4 + ^2 for 10 Rs ™ ~ r ~ 1 AU ( 20 ) 

5 D&M [r] 1.02 i?sun ^ r < 10 # S un , v 

5 m &nM 01 10 i? S un < r < 1 AU 1 J 



3.3 Global Electric Field 

If one considers the small mass of an electron compared to every other ion, then, one can 
derive the way demonstrated in Onel (2004), a global electric field E, shown in Figure 3, 
from the spherical, hydrodynamical Euler equation for an ideal, isothermal electron fluid, 
if only the electrostatic force density is taken into account: 



E[r] = — 

e dr 



In 



N[r] 



(22) 



Using Equation (22) one can calculate an electrostatic potential difference of 1.8 kV be- 
tween the photosphere (r = 1 -Rsun) and the Earth (r = 1 AU) by assuming the introduced 
density model Equation (18). 
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Figure 3: Left: The global magnetic field given by Equation (21) is shown. Right: The global 
electric held given by Equation (22) is illustrated. 

4 Discussions 

Within this Section two diagrams (Figure 4), one for the inner and one for the outer 
transport, will be explained. These diagrams are not related to any particular solar event. 
Their purpose is to show how those quantities are effected by global magnetic and electric 
fields, and local Coulomb collisions. Both diagrams in Figure 4 contain an altitude-time 
r[t], electron plasma frequency (from Equation (l))-time f[t], energy-time W[t] and a 
pitch-angle-time Q[t] := total M evolution for an electron. (Please keep in mind, that 
due to the chosen centre of mass system those quantities describe only the effective one- 
particle system electron-plasma background particle.) In addition, the energy amount 
belonging to the velocity component parallel (drift energy) and perpendicular (Larmor 
energy) to the magnetic flux density is plotted within the energy-time diagram. 

In both diagrams, based on usual our observations of type III burst origins (e.g. Figure 1), 
an electron at an initial altitude (of approximately r ~ 1.127 i?sun), which corresponds 
to an electron plasma frequency of 300 MHz, with a total initial energy amount of 50 keV, 
and an initial pitch-angle of 0t o tai,o — 0.25 i\ are assumed. This means, that the initial 
velocity components parallel and perpendicular to the magnetic flux density are equal, so 
that at the acceleration site it = to) the drift and the Larmor energies for the electron are 
exactly half of the total energy amount. 

As it can be seen from Figure 4 the electron moving toward the Sun transfers its drift 
energy into Larmor energy due to the magnetic (mirror) forces, while it loses energy due 
to Coulomb collisions. The energy gain caused by the electric forces is insufficient to 
accelerate the electron, which stops in its radial movement about 0.4 s later. The electron 
moving toward the interplanetary space behaves reversed: It transfers its Larmor energy 
into drift energy due to the magnetic mirror forces, while it also loses energy due to 
Coulomb collisions and the electric forces. Due to the decreasing particle density of the 
atmosphere and its high drift energy, this electron is able to leave the solar atmosphere, 
with less energy loss. Note, that the transport calculations are done 1 000 times for each 
transport direction. Afterwards all single results for each direction have been averaged 
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Figure 4: Two examples for the electron transport model introduced in Section 2 are shown. 
The diagrams are explained and discussed in Section 4. 

arithmetically in the shown Figure 4 to obtain a smooth plot. 

5 Summary 

A 1-dimensional model for the propagation of energetic electrons through a certain plasma 
background, with known average properties has been introduced. 

In case energetic electrons produced e.g. during solar flares can penetrate into the inter- 
planetary space, they can be measured in-situ, e.g. by the WIND spacecraft. In addition, 
energetic electrons are held responsible for the non-thermal X-ray radiation from the 
chromosphere, which can be observed by RHESSI spacecraft. 

Within this paper we performed a numerical study to characterise how the typical quanti- 
ties during electron propagation could be affected, by considering average models for the 
local plasma density, the global electric and the magnetic field. 

In our opinion, the main improvement of the deduced model, which easily can be modified 
to treat proton or ion propagation phenomena, compared with earlier common models of 
the community (e.g. Brown 1972; Bai 1982) is the innovative way of treating the pitch 
angle by introducing a binary dice. 
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